Network Synthesis Using Tchebycheff 
Polynomial Seriest 

By SIDNEY DARLINGTON 

(Maouscript received April 17, 1952) 

A general method is developed for finding functions of frequency which 
approximnte assigned gain or phase characteristics, within the special class 

of functions which can he realised exactly as the gain or phase of finite 
networks of linear lumped elements. The method is based upon manipula- 
tions of two Tchehycheff polynomial series, one of which represents the 
assigned characteristic, and the other the approximating network function. 
The wide range of applicability is illustrated with a number of examples. 

1. INTRODUCTION 

Network synthesis is tlie opposite of network analysis — namely, the 
design of a network to have assigned characteristics, as opposed to the 
evahiation of the characteristics of an assigned network. In general, 
there are spccKications on the internal constitution of the network, as 
weU as requirements relating to its external performance. A common 
form of the general problem is the design of a finite network of hnear 
himped elements, to jiroduce an assigned gain or phase characteristic 
over a prescrilied iuter\-al of useful freciuencies. The present paper re- 
lates to this particular form. 

In general, the restrictions on the network are such that the assigned 
performance cannot lie matched exactly. This gives rise to an appioxi- 
mation or interpolation pi'oblem. Foi' present purposes, the i>roblem i.s: 
to choose a function of fre([ucncy which matches the assigned gain or 
phase to a satisfactoiy accuracy, from that special class of functions 
which can \w i-ealized exactly with physical finite networks of linear 
lumped elements. The function of frequency may be defined in terms 
of network singularities (natural modes and infinite loss points). The 



t Proscntpd orally, in Ijcicfer form, ;iL the 1951 Western Convention of the 
Institute of Radio Jllnginccrs, and at the Symposium on Modern Network Sjjn- 
thesis sponsored by the Polytechnic Institute of Brooldyn and T!ic Ollico of 
N.ival ResearcSi, New Yorlv City, A]iril, 1952. 
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interpolation, problem may then be regarded as solved when a suitable 
set of network singularities has been obtained; for quite different tech- 
niques are used to design actual networks with these singularities. 

The interpolation problem may be attacked in a number of different 
ways; and a variety of different techniques are, in fact, needed to cover 
the wide diversity of practical applications. The present topic is a fairly 
general way of attacking the problem, based upon manipulations of two 
series of Tchebycheff polynomials. The two series represent expansions 
of two functions of frequency — ^one, the ideal assigned gain or phase, 
the other, the network approximation to the ideal. The interpolation 
problem may be solved in this way because it is feasible, as will be 
shown, to deteiTQine network singularities from arbitrarily assigned 
values of coefficients in the corresponding Tchebycheff polynomial series. 

The techniques to be described were deri\'ed originally from studies 
of the so-called potential analogy; but they can now be developed most 
easily %vithout reference thereto.! In a sense they may be regarded as 
extensions of familiar filter theory, using Tchebycheff polynomials, to 
more general gain and phase functions. The extensions, however, depend 
upon a number of new principles. Extensions of the filter theory apphed 
to more general problems have been noted in pubUshed papers; but 
those noted have not used the specific general approach employed here. J 
The wide applicabiUty of this general approach will be illustrated by 
specific examples. 

2. NETWORK AND TRANSMISSION FUNCTION 

It will be sufficient for our present purposes to limit the discussion to 
the general 4-pole shown in Fig. 1. The 4-pole may be active or passive, 
but it must be a stable finite network of linear lumped elements. E and 
V are steady state ac voltages, E the dri\dng voltage and V the response. 
The gain a and phase are here defined as the real and imaginary parts 
of log V/E, 

For a finite network of lumped elements, a -\- i(3 always has the fol- 
lowing form : 

a + i,3 = log 7C /" - ^i^l?' - P'i - • • (1) 



t TchebvchelT polynomials are related to potential analogue charges on el- 
lipses, as described in the author's paper "The Potential Analogue Method of 
Network Synthesis"^ Section 15. 

t For the most part, they have used the potential analogy, in such a way that 
TchebyclielT polynomials do not appear at all in general applications. For ex- 
amples, see methods of Matthaei^ Bashkow^ and Kuh^. 
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The "frequency variable" p represents, of course, ioj. The zeros p, of 
the rational fraction arc those values of p at which there is infinite loss. 
The poles p, are the so-called natural modes, or values of p at which 
response V can exist in the absence of dri\ing voltage E. The scale 
factor K determines the average level of transmission. The zeros, poles^ 
and scale factor together determine the gain and phase completely. 

For a physical stable network, the zeros and poles must meet certain 
well known restrictions, which are commonly stated in terms of loca- 
tions in the complex plane for frequency variable p. Withui these re- 
strictions, the zeros and poles can be subject to arbitrary choice, say 
for purposes of net\vork synthesis. 
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a+ i./3 = LOG V/e 
Fig. 1 — A general 4-pole. 

The symmetries required by the physical restrictions permit a and ^ 
to be represented separately as follows:! 

2a = log K' ^P'^' - P^^P^ - V"^ ■■• 

(Pl - P){p2 - p) ■■■ . 

i28 = log 



(p'l + p) ■ ■ ■ (pi' - p) • ■ ■ 

These expressions hold at all real frequencies, but only at real fre- 
quencies. 

3. TCHEBYCHEFF POLYNOMIALS 

It is functions of these special types which we are to synthesize with 
the help of Tciiebycheff polynomials. INIore generally, TchebychefT 
polynomials come in various forms, and may be analyzed in various 
ways. For oiu- special purposes, however, they take somewhat special 
forms (a little different from textbook definitions); and they are best 
analyzed in quite special ways. J It mil be simplest to start vnth arbi- 
trary definitions, to be justified later on by demonstrations of usefulness. 

t Tlie ])!iaae equation omits a possible 1S0° phase reversal, which is trivial for 
present, purposes. 

t For discussions of TchebychefT polynomials from other viewpoints, see 
Courant and Plilliert'', and also a paper by Lanczos^ on trigonometric interpola- 
tion. 
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Actually, the definitions must vary with the nature of the useful fre- 
quency interval, For the present, however, it will be assumed that the 
useful interval extends from w = to o), ; or more precisely, from 
o) = — ojc to +wc (in accordance with the symmetries of gain and phase 
functions). Useful intervals which do not include oj = require changes 
in the definitions, which \\'ill be taken up in Section 28. 

For our present purposes, Tchebycheff polynomials Tk may be defined 
as follows: 



p = to) = loic sni <t> 
Tk = cos k4>, k even 

Tk = i sin ^■<^, k odd 



(3) 



The first equation defines an auxiliary angle variable, 4>, in terms of 
which Tk is especially simple. The imaginary scale factor i, associated 
mth polynomials of odd order, simplifies later analysis. In addition, it 
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Fig. 2 — Tchebycheff polynomiala. 
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is especially appropriate for general network applications, because the 
odd ordered polynomials conliil)ute to the imaginary parts of complex 
network functions — such as i/3 in a + i&.f 

It is apparent from (3) tliat the TchebychefT polynomials become 
simply Fourier harmonics, if they are plotted against a distorted fre- 
quency scale — that is, against <^. This means that they must be ortho- 
gonal, over that particular range of frequencies which corresponds to 
real values of <^. From tlic relation between <^ and w, it is clear that real 
values of ^ cover the frequency interval between — Wc and +We , which 
is our useful interval. In other words, the interval of orthogonality coin- 
cides with the useful frequency interval. The corresponding interval of 
p is of course p = — lojc to -j-iuc ■ 

If a given function is plotted against ^, instead of w, it may be ex- 
panded in a Fourier series. Each term in the series may be replaced by 
a Tchebycheff polynomial, to obtain au expansion of a given function 
in terms of polynomials, for the specific useful interval w = — Wf to 
-\-Wc . Established techniques are available for expanding experimental, 
or other numerical data, in Fourier series, as well as actual analytic 
functions. 

In Fig. 2, some of the Tchebycheff polynomials are plotted against oj. 
The frequencies — Wc and +Wn are also indicated. Frequencies between 
these limits correspond to real values of the angle variable 0. If this 
part of the frequenc.v scale is stretched, in the proper non-uniform way, 
the \'arious "loops" not only \v.\ve the .same maximimi values, but also 
the same shapes. In other words, they become periodic. ^lore specifically, 
a stretch which changes the frequency scale into a <f> scale changes the 
plots into sin A'l^ or cos /,-<^. 

4. TRANSFORMATION OF VARIABLE 

An alternate to (3) may he obtained by relating a new variable, z, to 
4> by 

(4) 

Substituting z in the exponential equivalent of sin <t>, in the first equa- 
tion itf (3), gives an alternative definition of z directly in terms of p, 
namely: 

P = |(.-i) (5) 



t A smiill change in the definition of <ii would bring the definitions closer to 
convention, by replacing both sines by cosines (without altering Tk as a function 
of p). Tliis however, would complicate our later analysis. 
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Substitutions in the exponential equivalents of the other sine and cosine 
in (3) give: 



1 



r^ = ^ U + -J > '^ even 



1 



r. = ^ U' - rJ . ^^ odcl 



(6) 



Network applications depend upon the nature of the relationship be- 
tween the variable p, and the variable z. The relationship is illustrated 
in Fig. 3, which indicates corresponding contours in the p and z planes. 

Since angle ({> is real in the useful interval, z, as defined by (4), has 
unit magnitude. In equivalent coiiformal mapping terms, the unit circle 
in the z plane maps onto a segment of the axis of real frequencies in 
the p plane — namely the segment extending from p = —ic^c to +^coc . 
Hereafter, we shall say merely that the useful interval in the z plane is 
the unit circle, or | s 1 = 1. The real frefiuency intervals outside the 
useful interval map onto the imaginary axis in the s-plane. 

z-plane circles with radii other than unity map onto p-plane ellipses, all 
with foci at p = d=iuc . This is reminiscent of filter theory using Tcheby- 
cheff polynomials, and in fact such a filter may be obtained by spacing 
2-plane mappings of natural modes uniformly around such a circle. f 
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Fig. 3 — The complex planea for p and z. 



t The filter theory is developed in detail in a monograph by Wheeler'', which 
also includes an extensive bibliography. 
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5. Z-PLANE MAPPINGS OF NETWORK SINGULARITIES 

3-plane mappings of network singularities are also an essential part 
of synthesis applications. The mapping z^ of a typical zero or pole p^ is 
illustrated in Fig. 4. From (5), the analytic relation must be: 



V" = 



1 



(7) 



By its quadratic nature, there must be exactly two v^alues of z^ , corre- 
sponding to one Pa. The relation is such that replacing 2„ by —l/z„ 
leaves p, unchanged; and hence the two values of s„ must be negative 
reciprocals, each of the other. Thus, one mapping of p„ falls outside the 
unit 2-plane circle, and the other inside. 

A unique definition of z„ may be obtained by requiring that z„ must 
lie the mapping outside the unit circle. Then \ z„ \ > 1 hy dehnition, and 
■the complete definition of z^ may be: 



(8) 



I 2 J > 1 



This definition is unique provided network singularities p„ are excluded 
from that very special line segment of the real frequency axis which 
corresponds to the useful frequency interval, — Wc < w < +Wc (where 
] 2ff I would be exactly 1). 

We are going to solve the interpolation problem by choosing the z^ 
first, instead of the p-plane singularities p„ , after formulating the inter- 
polation problem in suitable s-plane terms. For this, however, we must 
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know what further conditions must be imposed upon the z^ , so that 
the corresponding p„ will meet the special conditions necessary for 
physical networks. A simple analysis of the definition (8) of z^ , and of 
the well known restrictions on the p, , leads to the following assertion; 

The 'physical restrictions on z„ are 
exactly the same as those oyi p„ . 

It is obvious, for example, that conjugate complex z, are necessary 
for conjugate complex p„ . Also, because \za\ > l,z^ dominates — l/Zc. . 
Tlien the sign of Re p., is the same as that of the Re z, , and p„ with 
negative real parts require z„ mth negative real parts, and so on. 

Thus the direct choice of z, is restricted in exactly the same way as 
the choice of p, , except for the additional general requirement | 3„ | > 1. 
The last condition imposes no important restriction on the corresponding 
p„ . Initially, it was adopted to make z, unique for any p„ (not at a 
useful real frequency); but this condition does also play an essential 
role in the 2-plane formulation of the interpolation problem. 

6. NETWORK GAIN AND PHASE IX TERMS OF Z 

A first step in the z-plane formulation of the interpolation problem is 
the formulation of the network gain and phase functions, (1) and (2), in 
terms of z. This is most usefully examined as a transformation of func- 
tional form, rather than as a conformal mapping. 

The gain and phase function (1) transforms as follows: The analytic 
relation between p and z is regular in the neighborhood of the singular- 
ities p„ of the network function. Therefore, there will be similar singu- 
larities of the transformed function at the 2-plane mappings of po , 
which are z„ and - \lz, . These singularities, and also suitable behavior 
at infinity, are exhibited by the following expression for a + i0 as a 
function of z. 

XI is used here to designate a product of factors of the type following it.f 

t The expression is readily confirmed in the following very elementary man- 
ner: For every factor M - |-V in (9), there is also a factor (l + — j . The 
product of the two TCiay be expanded as follows: 



( 



2„2 



(10) 
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If we define a new scale factor K^ by K^ — Kz , we may write (9) as 
follows : 



a -\- ip = log < 



L'^H)! 



n 1-4, 



/Same Rational 

\ Function in — l/z 



(13) 



Similar expressions for the separate gain and phase functions may be 
derived from (2) : 



[Same Rational 



2« = »«H^'jj/\_^\ \ Function in -lA 



1- ^ 1+4 

..,,, , I TT ^ TT ^° I /Same Rational 

,2^ = log m — r n — 7}i i.„,rti„„ i, _ 1/, 



z 

1+ -7 

2» 



2 
1 - ^ 

2ff 



(14) 



Equation (13) holds at all values of p and z, wliile (14) holds at all 
real frequencies. Simplifications of (l-i) should be noted, good for the 
useful interval only. When \z\ = 1, I/2 = 2*. Recalling also that log 
j X ^ is 2 log I .t; j, and similar elementarj^ relations, one obtains from (14) : 

When I 2 I = 1, 






a = log 



1-1, 1 + 4, 

)3 = Phasen- ill 7 

1 + -^ 1 - -7, 

2ff 2« 



(15) 



Comparison witli (ol and (7) gives: 



(' - 7) {' + £) = -„T. '" - ''■ 



Thus (9) ia equivaleut to (1) provided 



n - 



k: = K 



n - 



(11) 



(12) 
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7. THE POWER SERIES IN Z 

Our applications to network synthesis depend upon a correspondence 
which may be shown to exist between certain functions of z and certain 
power series in z. The functions of z may be formulated in terms of 
network singularities. The power series in z may be derived from the 
Tchebycheff polynomial series in -p representing the corresponding gain 
and phase. 

The Tchebycheff polynomial expansion of a gain and phase function 

may be written; 

a + z-^ = E G^'Fu (16) 

If a + i/S corresponds to a finite network, it may be represented by the 
function of z in (13). At the same time, Tk may be represented by the 
function of z in (6). With these changes, (16) becomes: 



log 



(,n(.-5)l 



n(i-i' 



/Same Rational 

\ Function in — Ifz 



(17) 



The logarithm of the product of the two rational functions, in z and 
— 1/s respectively, may be written as the sum of two logarithms. The 
series in sums of ^ and (-I/2)'' may be written as the sum of two series. 
Then 




, /Same Rational \ 
"^ ^^\ Function in ~l/z/ 



(18) 



=E^^/+E|(^y 



The above expression equates the sum of two similar functions, in z 
and - \jz respectively, to the sum of two power series, also respectively 
in z and — \jz. The theorem on which the synthesis methods are based 
asserts that the functions and power series in 2 and — Xjz may be equated 
separately, throughout the useful interval. That is: 



NETWORK SYNTHESIS 



623 



When \z\ = 1, 



log^ 


(,nO- 


-D' 


I'nO- 


-0 



log 



Same Rational 
Fuuction m —1/z 



= \ E cj' 



\i:gA— 



(19) 



The relation (18) does not, by itself, require (19) to be tiiie. (19) fol- 
lows from (18) if and only if tlie function of z has a power series expan- 
sion involving only positive powers of z, and the function in — 1/z has 
a power series expansion in — 1/z, mth the same coefficients. This added 
condition, howc\'er, is readily established for the useful interval. f 

Combining (19) and (16) yields a most useful relationship connecting 
the 2-plane mappings z^ , of the network singularities p, , and the 
coefficients Ck , of the Tchebycheff polynomial expansion of a + iQ: 



i: ^Ckz' = log /C 



n 1-5 



(20) 



n(i-Z' 



In more qualitative terms: 



The transformation from variable p to variable z 
converts an expansion in Tchebychef polynomials 
in p into an expansion in a power series in z. 

Thus, l)y working mth the z„ , in place of the p„ , one may use a 
power series sort of analysis in calculating, or in choosing, the coeflfi- 
cients Ck in the Tchebycheff polynomial series. 

The relations (20) refer to the combined gain and phase function. 
Exactly similar relations can readily be obtained, however, for gain and 

t A3 denned in (S), | z, 1 > 1. In the useful interval, \ s] =1. Hence 1 z/z„ \ < 1. 
It follows that log (1 — z/z,) has a power (MacClaurcn) series expansion in posi- 
tive powders of r, convei'sent on and within the circle | z | = 1, Finally the firat 
logarithm in (19) may be expressed as a sum of logarithms of this simple type, 
each of which may be exjianded sejiarately. Substituting —l/z for s maps the 
unit circle onto itself. It follows that the second logarithm in (19) has an expan'- 
sion in positive powers of —1/z, in the useful interval, provided the first loga- 
rithm has an expansion in positive powers of z; and the coefficients in the two 
series will be the same. 
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phase separately. These may be derived from (14), and take the form: 



Y.Cuz" = logAl 



i& = E CuTu 



z- 



mi-u 



n(i-|= 



1- - 1+^7 

Ec.3^ = iogn^^n — f 

1 + - I - -n 

Z„ Za 



' k, even 



(21) 



• A:, odd 



(The absence of factors | in E '^'a^*- as compared with (20), reflects the 
factors 2 associated \dth a and ,3 in (14).) 



8. REPRESENTATION OF ASSIGNED GAIN AND PHASE 

In synthesis problems, the networlv gain or phase, a or /3, is to ap- 
proximate an assigned (ideal) gain or phase, say a or ^. To make effec- 
tive use of the 2-plane analysis, in network synthesis, we need to describe 
a and ^ by relations analogous to (20) and (21), which express a and Q 
in 2-plane terms. These relations, while similar to (20) and (21), must 
take a more general form (since a or ^ need only be approximately the 
gain or phase of a finite network). For our present purposes, the ap- 
propriate relations are those noted below. 

Let a. + ?^ be any function of p which has the following properties: 
It must be analytic throughout the useful interval. Further, there are 
to be no singularities within a (p-plane) distance e of the useful interval, 
where e is finite (but may be small). Finally, at real frequencies, a and 
i$ are to equal respectively the even and odd parts of a -f- ij3. 

Under the conditions stated, a + ifi may always be expanded in 
terms of our Tchebycheff polynomials Tt . Let E ^^Tk be the expansion. 
To obtain a parallel to (20), we may form_(arbitrarih') a power senes 
E iCkz''. Then we may define a function Ii(z) by identifying log R(z) 
with the power series. All this adds up to the following, comparable to 
(20): 

E hC^z' = log R(z) 
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The functions of z have the following properties: Because of the mild 
restrictions, which we have imposed on the singularities of a + i^, the 
series zl Ci-z'' dehnes a function which is analytic within, and on the 
circle \z\ = 1. Then R{z), also, is analytic within, and on the circle. 
Further, ii{z) has no zeros anywhere in the same region. {R{z), however, 
may be more general than the rational fraction in (20).) FinaUy, because 
of the even and odd .symmetries, required of a and i^, (22) may be 
broken into the following parallels of the equations (21): 



EcJ- = \og[R{z)R(-z)] 
R(z) 



E Cj = log 



R(-z) 



k, even 



/i, odd 



(23) 



In some applications, it is po.ssible to express R(z) in closed form. In 
all applications, it is possii)le to expand R(z) as a power series, convergent 
in the region \ z \ ^ 1. The same is true of 1/R(z), since there are no 
zeros in the region. Coefficients of cither series (R(z) or 1/R{z)) may 
readily he calculated by means which we shall e.\-amine a little later. For 
the present we shall .say merely tirat R{z) is a known function, corre- 
sponding to an assigned a -\~ i^. 



9. A DESIGN- CRITKRIOX 

When the gain and phase fuiiclion, a + ii3, is to ajiproximate a + 7'^, 
the error in the approximation is (« — «) + i(i3 — $). The error may 

be expressed in terms of z by taking the difference of corresponding 
equations in (20), (22). The difference of the logarithms may be ex- 
pressed as a single logarithm of a ratio. Alternatively, and also more 
conveniently for our later purposes, it may be expressed as the negative 
of the logarithm of the reciprocal ratio. When this is done, 



(a - a) + 7(|3 - ~&) = E (Cl. - C,)n 



EKc. -c,)/= -iog{ ~ 



U i-^ 



n 1- 



■Riz) 



(24) 



Consider the following arbitrary requirement, as a design criterion: 
The series E (-^Tk is to match exactly the .series E CkT,, , through 
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terms of order m. If both seiies have converged to small remainders 
when k = m, this criterion will surely make a + t/3 a good approxima- 
tion to a -\- ip.\ In terms of the coefficients, the criterion requires: 

C, = C, , k ^ m (25) 

If (25) is applied to the second equation of (24), the power series is 
zero through terms of order m. In other words, the logarithm, equated 
to the series, will approximate zero in the power series, or "maximally 
flat" manner, to order m. The logarithm is zero when the expression in 
brackets is unity. Further, the logarithm will approximate zero in the 
maximally flat manner when, and only when the bracket approximates 
unity in the maximally flat manner. Thus a condition which is equivalent 
to (25) is the following: 



J- —A !l/ .R{z) = 1 + e„+i z-"" + 6^2 z-"-' ■ ■ ' (26) 

This may be represented symbolically by 

in(i- 



z 

Z. 



^-n(i-5 



z« 



H{z) = 1 (27) 



where = is used to indicate equality through power series terms of 

order m. 

^Vhen (27) is applied to network synthesis, the singularities z„ , and 
scale factor K^ are the unkno-\vns, while R{z) is known. If vi is equal to 
the total number of z, , (27) will determine the network function com- 
pletely. When m is smaller, (27) will furnish m + 1 relations between 
the network parameters (including the zero oixler condition), which 
may be combined with specifications of other sorts. Since (27) is equiv- 
alent to (25), this procedure amounts to the determination of network 
parameters which will >4eld assigned values of the coefficients, Ct = Ck , 
h ^ m, in the Tchcbycheff polynomial expansion of « -|- i&. 

Equation (27) applies when both gain and phase are to be approx- 
imated. For approximation to gain only, or to phase only, similar rela- 
tions may be derived from (21) and (23). Only even ordered Tchebycheff 

t When both residues are relatively large, the approximation may_ still be 
good, for the remainders may be quite similar, and the error will be their diner- 
ence.' In practical ap])lications, thi.s is a not unconimnn situation. 
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polynomials contribute to gain. The following condition turns out to be 
the equivalent of C2k = C2k , /^ ^ m: 



i_n(i- 



z' 



z'J' 



""^nd- 



■R{z)R{~z) = 1 (28) 



where = means approximation in accordance with a power series of 
even ordered terms, througli order 2m. Correspondingly, only odd 
ordered Tchebycheff polynomials contribute to phase. The following 
condition is equivalent to Cn—i = Cik-i ,h = 1 to m: 

z z 

L -\- II i — / 

Za Za 

The remaining sections (except the last) develop in more detail the 
application of 2:-plane techniques to more specific synthesis problems, of 
various sorts. Most of these (but not quite all) are based directly on 
(27), (28), or (29). The exceptions use a modification of (28), in which 
the function of z on the left is retained, but mth the zeros and poles 

me 

adjusted for a different kind of approximation to unity, = but not = . 
In all cases, miity is approximated with one of the functions appear- 
ing in (27), (28), (29). It will be convenient to use H{z) to represent the 
error in the approximation, or departure from unity. When gain only is 
of interest, the function in (28) is used, and H{z) is defined by: 



i n(i-f-. 
^^n(.-4 



'U{z)li{-z) = 1 + H{z) (30) 



In developing the specific techniques, we shall start vAi\\ a very 
definite, rather special examiile, in order to illustrate the techniques 
witli specific operations. This will be discussed in considerable detail in 
Sections 10 through 14. Thereafter we shall examine how these specific 
operations may be generalized, in a number of different respects. 

10. AN INTRODUCTORY EXAMPLE 

The example which has been chosen for detailed discussion is the 
equalization of the gain distortion produced by two resistance-capacity 
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type cut-offs. The equalization is to be accomplished with a network 
which has n natural modes, but no finite frequencies of infinite loss. 
(This is simply one of the arbitrary specifications which define this 
problem.) 

The two cut-offs may be due to circuits or devices at two different 
points in a communication system, which may be represented schemat- 
ically as in Fig. 5. Their effect can be described in terms of tw-o assigned 
natural modes. Two assigned modes are assumed, instead of only one, 
because a single mode would make the problem too shnple. Our present 
purposes will be served well, however, if we simplify the problem by 
requiring the two assigned modes to be identical, say at p = pa ■ 




rig. 5^A system with two resistauce-capacity type cut-offs. 

The two natural modes would be cancelled completely by two infinite 

loss points at the same location in the p plane. A network with two 
infinite loss points, however, is not physically possible unless it has also 
at least two natural modes; and the natural modes will have to introduce 
distortion of their own. Thus no finite network will give perfect equali- 
zation of unwanted natural modes. Sometimes it is desirable, in practice, 
to use an equalizer configuration which produces no finite frequencies of 
infinite loss, the entire equalization being accomphshed by a suitable 
choice of its n modes. Configurations of tliis sort are illustrated in Fig. 
6. Thus, our simple illustrative problem, though chosen to introduce 
principles, is also of some practical interest. 

The exclusion of finite frequencies of infinite loss simplifies the repre- 




rUT^ 



<innr^ 



^lyw^ 




Yig, 6— Coiifiguratioua which produce no iiuitc frequencies of infinite loss. 
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sentation of the network gain a. In {21 ), tlio z^ rorr('s])i)ii(l to finite fre- 
(]neiu'ies of infinite loa«, and are to be omitted wlien thei'e are to be 
iiatnriil modes only. What i.s left is the logarithm of the reciprocal of a 
polymimial, which is of eonrsp the negative of the logarithm of the 
p()lyn{)inial itself. Thus a may lie described as follows, for this particular 
application: 

E CW = -lof^ K: n f 1 - -2") ^ = 1, • ■ -, n. (31) 



(For convenience, z^ has boon written for z^ , and K^ has been redefined 
to avoid the l/Id, required if it is defined as in (21).) 

The assigned gain a is even more special. In this particular problem, 
a has most of the properties of a network gain a. Specilically, it is the 
negative of the gain to be e(|ualized, wliich in fact corresponds to a 
finite network. As a result, R(z) of (22) may be expressed in closed 
(rational) form. (Later on, we shall modify the methods appropriate for 
this very special situation, so that R{z) need l)e representable only by 
series.) 

The specific representation of our present a may he very similar to 
the representation of a in (31), as follows: 

a = z^ C-2kT<ik 
Ec^2^'' = +logX;(l -0 ^^-^ 

(Both (31) and (32) apply only to the useful interval, i ^ | - I.) 

The constant Zn is the z-plaiie mapping of the assigned unwanted 
natural modes oX, p = p» , aiid may Ije calculated tlierefrom liy (8). In 
(32), j(i detci'mines the C--k- , which in turn determine a. The constants 
Zb , in (31), are the 2-planc mappings of the arlntrary natural modes of 
the e(iualizer. They are to hv adjusted to make a approximate a. Then 
the network natural modes /;„ may be calculated from them, by means 
of (8). 

TakiTig the difference of corresponding equations in (31) and (32) 
gives the following, analogous to (24): 

a — a = ^ [C-ik — C'-k)T2k 

E ic,^ - c..)z'' = -log [k^k: (1 - ly n (1 - I)} ^^^^ 

This differs from (24) in two regards. It relates to the gain error, {a — a), 
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without regard to phusp. It reflects the more specific functional forms of 
(nu" present a and a. 

The formulas show that the coefficients in the Tchebycheff poly- 
nomial expansion of our present a - a are fixed by the logarithm of a 
polynomial in z', of degree n + 2. Since the Tchebycheff polynomial 
series is simply one representation of the function a — a, this means 
that a - a itself is determined by the polynomial in si Out of the n + 2 
zeros, in terms of z, n are subject to arbitrary choice, but the other two 
are required to be at z' = Zo . 

To arrive at a useful choice of the zeros, one may start with the ex- 
panded form of the polynomial, which replaces the second equation of 
(33) by: 

E (C,>^ - C,,)^'" = -log \Iu + fC^z' + ■ ■ ■ K.+2z"'^'\ (34) 

All but two of the coefficients K^ may be assigned arbitrary values, pro- 
vided the remaining two are then adjusted to give the reciuired two 
zeros at z' ^ zl . The corresponding zeros zl may then be found by 
ordinary root extraction methods. 

The coefficients may lie chosen in such a way that the complex poly- 
nomial approximates unity, when | ^ j = 1. Then the logarithm approx- 
imates zero, tiie coefficients in the power series (34) are small, and since 
these are also the coefficients in the Tchebycheff polynomial series in 
(33), a — a is small. 

II. TCHEBYCHEFF POLYNOMIAL SERIES MATCHED THROUGH tl TERMS 

A special choice ()f coefficients, which meets these requirements fairly 
well, is the choice determined by (28), with m = n. The function on the 
left side of (28) is here the polynomial in (34). For our present purposes, 
therefore, (28) becomes: 

{TC-o + AV + • - ■ lu^-^z'-'^'] = 1 (35) 

This requires A'o = 1, and Kk = for k = lion. Then i?„+i and i?„+2 
must be adjusted to give the two required zeros at z' = zq . This gives: 

Ao + • . ■ K^^^'"^' = 1 - ( » + 2) (iy"""' +(«+!) (jj^^ (36) 

In accordance with Section 9, this special choice of coefficients corre- 
sponds to a match of Tchebycheff polynomial series, a to a, through 
terms of order 2n: 

C,, = Cn. /v ^ ^ (37) 
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Fig. 7 — Error when four natural modes equalize two natural modea at 
yn = — 0.75wr. 

The a(itual accuracy may be calculated from the final zeros and poles, 
from non-zero terms in the expansion of o: — a, or by an analysis in 
terms of z which will be described later. 

A sample plot of a — a is shown in Fig. 7. This corresponds to an 
equalizer of four natural modes, compensating for an initial loss which 
rises to about 8 db at the top useful frequency, or a distortion of ±4 db 
about the median loss. Residual errors are order of ±0.06 db. 

A little later we shall return to the question of accuracy, to take up 
methods of estimating what can be done with other numbers of arbi- 
trary natural modes, and other values of the assigned modes. First, 
however, we should investigate whether the network singularities z^ 
determined by (36) meet the other necessary conditions. 



12. PROPERTIES OF %„ 

In the first place, | 2, \ must be >1. Otherwise, the function of z in 
(31) will have no power series expansion over the useful interval \z\ = 1 ; 
and (31) will not, in fact, determine the gain a over the useful interval. 
It turns out, however, that the condition does not give trouble in the 
synthesis of natural modes, when there are no arbitrary freciuencies of 
infinite loss. This may be demonstrated by the argument outlined below. 

The Za are zeros of the polynomial in (34), which we have given the 
special form (36), by applying (35). A function theoretic test for I s^ 1 < 1 
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makes use of the contour in the complex plane for the polynomial, corre- 
sponding to the 2-plane circle \z\ = 1. (This is like a Nyquist diagram 
except that the contour for the variable, z, is different.) There will be 
\z,\ < 1 if and only if the contour for the polynomial encloses the 
origin. 

Now the polynomial in (34), and (35), is merely a special case of the 
function on the left in (28), and (30). For this special case (30) becomes 

Y.K,z'' = 1+H(z) (38) 

The polynomial cannot enclose the origin without passing through some 
negative real value. But this requires an | H{z) \ > 1, at some point on 
the contour in question, \z\ = 1, which happens to be also our useful 
interval. On the other hand, a — a = Q when ^ ^kZ^'' = 1, and Hiz) 
is in the nature of a correction term, which is small in the useful interval 
when a — a is small. 

The conclusion is: There will be no j 2, | < 1 unless the approxima- 
tion, a to a, is so poor that a — a exceeds several db in the useful inter- 
val. 

Besides the requirement \z^\ > 1, the z, must meet physical restric- 
tions, which we found to be the same as those limiting the natural 
modes p, . The 2, may be calculated as follows: The zl are roots of the 
polynomial in (35), in terms of 2^ All the roots in terms of z^ are zl , 
except the two required roots at zl , which correspond to assigned gain a. 
Each z„ is a square root of a 2^ . There are two possible square roots, 
however, differing only as to sign. As far as gain a is concerned, either 
choice of sign is permissible; for a depends only on zl . For a physical 
network, however, the choice must be such that Re Za < 0, This choice 
is possible if, and only if V^ has a finite real part. A pure imaginary 2„ 
corresponds to a negative real zl , and thus negative real roots in terms 
of 2^ are excluded by physical considerations. 

Table I lists both zl and 2„ for a number of values of n. When n is 
even, all roots are physical. On the other hand, when n is odd, one root 
is always non-physical. In a sense, an odd n is not really appropriate 
for the present illustrative problem, ^\'ith any physical design. An odd n 
must necessarily bring in a real natural mode, which merely increases 
the sort of distortion we are trjnng to equalize — that is the distortion 
due to unwanted real modes. 

The follomng argument substantiates the suggestion, and also illus- 
trates manipulations of a sort which are frequently useful in more gen- 
eral applications: The highest order coefficient in (34), i?„+2 , may be 
set aside for adjustments to satisfy physical requirements. The rest of 
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Table I — Z-Planc Natural Modes jor Equalization of Two Identical 

Unwanted Modes 

1 -.5000 ± i .7071 Non Physical 

2 -.3333 ± i .4714 ±(.3492 ± i .6747) -.3492 ± i .6747 

3 -.6059 ± i .7784 Non Phyaieal 
-.0720 ± ( .6384 zb(.5340 ± i .5977) -.5340 ± i .5977 

4 +.1378 ± i .6782 ±(.6441 i i" .5264) -.6441 ± i .5264 
-.5378 ± i .3582 ±(.2328 ± i .7695) -.2328 ± i .7695 

5 -.6703 ± i .8187 Non Physical 

+ .2942 ± i .6684 ±(.7157 ± / .4670) -.7157 ± i .4670 

-.3757 ± i .5701 ±(.3918 ± ( .7275) -.3918 ± i .7275 

the coefficients may then be choseu to eliminate terms from the series 
S (^at ~ Oik)T2k , representing a — a, subject to the previous condi- 

nt {n— l)e 

tion that two zeros must bo z' = So . This replaces = by = , in 
(35), and changes (36) to: 

^"^ ^"^ , (39) 

+ A„+22 (2o - 3 ) 

If n is odd, all the roots z^ can be physical only if i?n+2 is negative. 
On the other hand, any finite negative R„+i leads to a larger error, 
a — a, than Kn+t = 0. Reducing K„+2 to zero is the same us reducing 
the degree of the polyiiomial by one, which amounts to reducing n by 1, 
from an odd to the next smaller even integer. In other words, a physical 
design mth an odd number of natural modes is less effective, for the 
present application, than a simpler network, with the next smaller even 
number of modes. 

Note that the z, in Table I are proportional to Zo . This means that 
root extraction methods need be used only once for each value of n, 
after which the roots may be quickly adjusted for any v^alue of Zq , 
corresponding to any assigned value of the two identical modes, po . 

13. ACCURACY 

The accuracy of a completed design can be checked by calculating a 
from the natural modes p„ , and comparing a with a. It is important, 
however, to have at least some information about accuracy in advance 
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of the detailed calculation of the p„ . Otherwise, it may be necessary to 
carry out several designs, iu all detail, in order to obtain one satisfactory 
design. 

The needed information about accuracy can in fact be obtained from 
the error function H(z), which we formulated for general gain applica- 
tions in (30), and for the present application in {38). The analysis 
which yields (15) may be used to obtain a very similar expression for 
a — a, in which R{z)R{ — z) appears in combination with the rational 
function of z from (15). It may be expressed in terms of the error func- 
tion Ii{z) of (30), as follows: 

a- a= -log I 1 + H{z) I (40) 

When H{z) is zero, a — a is zero. When H{z) is small, a — a depends 
on phase li{z) as much as on ] H{z) \ . When H{z) is a positive real, 
a - a is negative. When H{z) is imaginary, a - a is very small. When 
H{z) is a negative real, a — a'ls positive. When H{z) is complex, \a — ot\ 
is always smaller than with a real H{z) of the same magnitude. The last 
statement may be expressed as follows: 

- log Si + I H{z) |} ^ a - a g - log (1 - 1 H{z) \\ (41) 

The left hand relation is an equality when phase H{z) is an even num- 
ber of T radians; the right hand side, when it is an odd number of ir 
radians. 

In the useful interval, where z = e'*, the H{z) corresponding to (36) 
is as follows: 



n + 2 

mz)\ = -i^;^ 



Zq 



1 _ J^, ,.* I (42) 

{n + 2)2,, I 

f " + 1 2.-*\ 

phase H{z) = tt + (2/i -^ 2).^ + phase |l - ^^ ^ ^^.2 e | 

As to varies from to w^ , <^ varies by - radians. The corresponding 

phase of ^(2) varies by {n -\- l)ir radians, which means that H{z) is 
successively positive real, imaginary, negative real, imaginary, through 
n -\- I half cycles. This accounts for the oscillatory nature of the a - 5 
curve, illustrated in Fig. 7. 

The amplitudes of the oscillations are fixed by | H{z) \ , which varies 
relatively slowly. Specifically, the two logarithms in (41) determine 



NETWORK SYNTHESIS 



635 



envelopes, between which the actual error curve oscillates. These are 
the dashed lines in Fig. 7. 

The maximum error, in the useful interval, is determined by the 
maximum vahie of the envelopes, i.e., 



2o 



1 + 



71+2 



2 id 



(43) 



This function is plotted in Fig. 8, for various values of n. The abscissae 
"distortion before equalization" represent distortion relative to the 
median loss in the useful interval, or one half the total variation in the 
interval. (This is a function of the top useful frequency Wc , relative to 
the assigned natural mode pa ; and (7) makes Wc/po a simple function of 
zq .) The figure is convenient for estimating the values of n needed for 
specific applications. 

The various ripples in a — a do not all have the same amplitude, (43). 
For some values of n and z^ , the amplitudes are almost uniform; for 
others they are (juite varialile. A measure of the variability in ripple 
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amplitude, across the useful interval, is: 

I H jz) U... _ {n + 2)zl + (n + 1) 
jH(2)Lin. («+ '2)zl - in-h 1) 

14. APPROXIMATION IN THE TCHEBYCHEFP SENSE 

The above analysis suggests a way of improving the design deter- 
mined by (37) (or the equivalent, (36)). An optimum a - a is commonly 
one which has the following properties, in the useful interval: 

A Tnaxi7nu7n number of '^ripples," 
all maxima of \ a — a\ equal. 

(This usually minimizes the largest departure in the useful interval, 
thereby jiekling an "approximation in the Tchebycheff Sense.") Since 
the variation in phase H{z) determines the number of ripples, while 
I H(z) I determines the amplitudes of the ripples, the above conditions 
will be met if H(z) has the follomng properties, in the useful interval: 

Phase H{z) as variable as possible, ,^* 

I H(z) 1 constant. 

These conditions may be regarded as alternative design criteria, re- 
placing Cik = Cih ■ They can in fact be applied to our special example, 
and also to certain other special problems which will be noted later. 
For more general applications, a suitable H(z) can be defined, but no 
reasonably simple procedure has yet been found for calculating the re- 
quired constants. (The difficulties will be particularized in a later 

section.) 

For the present example, (38) may be used to replace (34), and hence 
also the second equation of (33), by: 

E (C'2. - C,,)z'' = - log [1 + 7/(^)1 (46) 

(33) requires H(z) to be a polynomial in z', of degree n + 2, with two 
zeros of [1 + H{z)] at / = Zo . The object is to find an H{z) of this 
sort, which also satisfies (45), at least to a good approximation. 
The following H(z) does in fact exhibit the required properties: 

H{z) = Gz [i-J/z^] ^^^^ 

The function is a polynomial because the factor [1 — /""'"V^^""'"^ is 
divisible by (1 — J/z^)- The constants J and G are to be chosen to 
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give the required double zero of [1 + H(z)] at z' = zl . One value of J, 
so detei-miued, is i-eal and of order 1/zl . This is the appropriate solution. 
Then | ./"^V^'"^"' 1 i« of order 1/zl"'^* , when \z\ ^ 1. This suggests the 
following approximation in place of (47): 

HU)^Gi'--l^_ (48) 

The approximation is at least as good as 1/20""*^*, compared witli 
unity, both in the useful interval and in the neighborhood of the singu- 
larities Zo , and z„ . This means that the approximation can be used: in 
estimating the error a — a (in the useful interval), in calculating ./ and 
G, and in finding the roots z„ . 

In the useful interval, \z\ = 1, and therefore 1/z = z*. Then 
(1 — J/z') is (1 — Jz')*; and their ratio has magnitude unity. Thus 
I ^(2) ! = I G I , iu the useful interval, to order of 1/zl"^* compared 
with unityf. With \ J \ < 1, phase H{z) varies over the useful interval 
to the same extent as the phase of ^"""*'^J Fig. 9 illustrates the difference 
ill a — a, as determined by (42) and (48). These curves, however, are 
for single values of n and 20 ; and the improvement obtained by using 
(48) would be different with different values of n or Zo . 

The values of J and G, determined from (48), and the requirement 
that [1 + H(z)] must have two zeros at z = zl , turn out to be: 

n + 1 1 2 



J = 



n-\-2zl , , 1 



/(-«)• 



l + 7^+/l/(l-^f + 



G = - 



zo V \ nj {n + 2fzt (49) 

1 1 - J/zl 



20'"^" 1 — Jzl 
Note that this / is iu fact smaller tlian l/zl . 

15. GENERALIZATION 

The several sections preceding describe a quite specific example, as 
an introduction to synthesis applications. The next several sections de- 
scribe how the specific methods of the example may be generalized, in 
several respects. 

First, the ideal gain, a, is generalized, so that it need not even liave 
the sort of functional form associated \\-ith finite networks. Then, the 

t TIh; I H(s) I (iclcrmined by (42) is constant only to order l/zl . 
t Tliis is tlie most, we can oxpeL-l, when we have n singularities, which can pre- 
vent the dominance of only lower order terms, through s'". 
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Fig. 9 — CompariBon of design procedures — four natural modes equalizing two 
natural modes at pa — — t^ we. 

approximating network gain is generalized, by introducing arbitrary 
frequencies of infinite loss, in addition to the arbitrary natural modes. 
The methods are also modified for approximation to an assigned phase, 
instead of gain, or to phase and gain simultaneously. Finally, the anal- 
ysis is modified to permit useful intervals of the "high-pass" type, or 
(in the case of gain simulation) of the "band-pass" t3rpe. 

16. APPROXIMATION TO A GENERAL ASSIGNED GAIN a 

If we now permit the assigned gain a to be general, in the sense of 
Section 8, we must return to the formulation: 

ZC^kz'" = \og[R(z)R(-z)] 
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If we retain simulation with a iietwoik which has n natural modes, and 
no frequencies of infinite loss, we must retain the formulation: 

The corresponding formulation of the error is {in place of {33)): 

z (Cu - c,,)z''' = -log [ai n (i - 1)-^(^)^(-2)] 

For Cafc = C-ik ,k^m, the following special case of (28) is now required: 

I- -.j-Rizm-z) = 1 (53) 

Now the reciprocal of R{z)R{ — z) has a power series expansion, in the 
region of interest. (Recall Section 8.) 

It follows that {53) may be multiplied by this quantity, without 
damaging the equality of power series coefficients. In other words (53) 
is equivalent to : 

Let Kk be the coefficient of z''' in the polynomial expansion of the left 
hand side; and let Kk be the coefficient of z^ in the infinite series ex- 
pansion of the right hand side. Then, 



J" 






(55) 



R{z)R{-z) 
Substitution in (54) gives 

K, + K,z' + . ■ . Kj" 1' E Kkz'" (56) 

In other words, 

Kk = Kk, k ^m (57) 

These relations are directly applicable to network synthesis, provided 
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the coefficients Kk can be calculated. Formulae for their calculation 
may be derived in the following way: 

If (55) is substituted in (50), the result is:t 

a = Y^ C2,T2k (68) 

When 3 = 0, the second equation reduces to: 

Ko = e-^o (59) 

If the functions of z are differentiated, a simple rearrangement gives: 

j:kKk/'^' = [-Z^^,..^'"][E^^2'*l (60) 

The right hand side may be expanded as a single power series, and then 
like powers on the two sides may be equated separately. The result is: 

Ki = — CiKa 

2K2 = -CjCx - 2CJC0 (61) 

3X3 = -C2K2 - 2CiKi - 3C^o 

Synthesis calculations may now be carried out in the following stages. 
The assigned gain a is expanded as a Tchebycheff polynomial series, to 
determine coefficients Cst , say through order k = n. The equations ((il) 
are then used to calculate coefficients Kk , also through order n. Each 
successive coefficient is computed in terms of those previously deter- 
mined. Note that the Kh , k S n, are fixed by the same number of 
Cik — that is, orders k ^ n. 

Equation (57) is now applied to identify Kk with Kk , k ^_ m. If all 
the network degrees of freedom are to be used to get Cat = C2k , index 
m = n, and (57) determines the polynomial in (55) completely. Other- 
wise, m < n, and coefficients K^+i to K^ are to be adjusted in accord- 
ance with specifications of other kinds. When all the Kk have been 
determined, the singularities z^ are found by root extraction methods, 
applied to the right hand side of the first equation of (55). 

The previous example might have been carried out in these terms, 
but happened to be simpler in the terms used. If (32) is regarded as a 
special case of (50), and if (32) is simplified (for purposes illustration) 
by using K, = 1, the corresponding R{z)H{-z) becomes simply 

1 -^l.ThenE^-^-^' 



IS 

3o. 



t We may think of these equations as defining an infinite network, with na- 
tural modes only, which would match the assigned gain a exactly. 
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(62) 



thus Rk and Kk become 



r^ fc + 1 

K-k = ■ 12 r~ , fc = to CO 

/^ + 1 

Kk = ~ 'lik ' 1 A; = to ?i 

2q 



(63) 



If these Kk are used to evaUiate the polynomial on the left hand side of 
(53), in accordant^e with (55), and if the polynomial is then multiplied 
by the above special R{z)R{ — z), the result is exactly (36). 

The error function H{z), of (30) and (42), may now be defined as 
follows : 



^' n (l - -^ R{z)R(-z) = 1 + H(z) 



(64) 



The error a — a is again : 

a - a = - log I 1 + H(z) 1 (G5) 

If (64) is used to express (53) in terms of H(z), a "1" may be sub- 
tracted from each aide of the relation, to get 

H{z) = (66) 

When m = n, this requires an H(z) of the following form, in terms of 
the coefficients Kk derived from R{z)R(—z): 

^ 2n+2 , ^ 2n-H _j_ 

H{z) = -'^^ +%t^_^L^ (67) 

17. CHAEACTERISTICS OF Z, 

As in the previous example, | z„ ] > 1 when the approximation, a to 
a, is at all rcasonaljle. The Za arc again zeros of 1 + //(z); but now 
H(z) is defined by (64). If [ H{z) \ < 1, when \z\ = 1, there will be 
the same number of zeros of 1 + H(z) as poles, in the region \z\ < 1. 
Any poles would liave to be poles of R{z)R(—z). In Section 8, we noted 
that this function is regular in the region [ 2 | ^ 1. Hence, there will be 
no poles, and there will be no [ 2^ | < 1, under ordinary accuracy con- 
ditions. 

As before, the Za can bo chosen in accordance with the physical con- 
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ditions, provided none are pure imaginaries. Since an imaginary z„ is a 
negative real z'g : 

There must he no negative real z„ . 

There will be no negative real zl if the polynomial ^ K^z^'' in (55) is 
non-zero at all negative real z. If the requirement is violated, initially, 
one or more Kk , of the highest orders, must be modified. Graphical 
methods are likely to be useful for this, combining plots of the original 
polynomial, and proposed changes. An approximation of the form 
C^k = Cik ,k -^ m, will still be realized, but with m < n. 

The error function corresponding to m < n is as follows, in place of 
(07): 

H{z) = tA^+i- J:»^i)^ _--^--K»-^.^ (68) 

18. ACCUR.\CY 

The accuracy of match again may be estimated by the means of (41), 
using H{z) of (67) or (68). H{z), however, may not be so easily cal- 
culated as for the previous example. 

A simpler but less reliable estimate of accuracy is furnished by the 
error in the first unmatched coefficient in the Tchebycheff polynomial 
series. If lU = Kk through k = m, the leading terms in various series 
are as follows. First, from (64) and (68), 

X^ n (l - I) Ri.m{-z) = 1 + ^""~/""' z"*' ■ ■ ■ (69) 
using this in (52) gives: 

E iCn - C..)^" = -log [l + { "'•" l^"'' ) ^'""' ■ ■ ■] (7 
Then, from the properties of logarithms, 

E iCu - C.k)z'' = - ^^±^^^^ .-- . . . (71) 



0) 



Consequently (also from (52)): 






K 

;eri 
placed by — Tim+i . If m = n, the same equation holds with K,„+i = 0. 



This is the same as the leading term of H{z), except that z^"""^^ is re- 
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The coefficient '""'"^ "^ in (72) is a sort of average of the enve- 

lopes of the ripples in a — a. The variabiUty of the envelopes, across 
the useful interval, depends upon higher order coefficients, in compari- 
son with the leading term. Calculation of higher order coefficients is 
relatively complicated. 

19. APPROXIMATION IN THE TCHEBYCHEFF SENSE 

The criteria (45) carry over to general assigned gains, as conditions 
on H{z) which, if i-ealized, are usually sufficient to establish approxima- 
tion in the TchebychefT sense. For this purpose we must use the Ii{z) of 
(64), rather than (07) or (G8) (which correspond explicitly to Cst = C^k , 
k ^ m). In terms of the polynomial and scries representations of (55), 
the H(z) of (64) becomes: 

The following somewhat special problem is easily solved, in these 
terms, and has a direct beiiring on various quite different synthesis 
tcchni(|ues: A network is to be designed which combines the functions 
of an etiualizer or simulator, with those of a filter, or selective network. 
In the useful interval, an assigned gain variation a is to be approximated 
in the TchebychefT sense. At higher frequencies, there is to be a rapidly 
increasing loss, or "sharp filter cut-off." The number of natural modes, 
n, is to be more than sufficient to match a to the required accuracy, in 
the absence of a selectivity requirement, the latitude being used to pro- 
duce the required sharp cut-off. In particular, n is to be large enough so 
that an n term match of Tcheljycheff coefficients produces errors that 
are negligilile compared with those accepted as a price of the sharp 
cut-otY. 

On the above assumption of an ample n, the infinite series ^ Ki.-z''' 
in (73) may l)e truncated after the term of order n, and the errors due to 
the truncation may be neglected iyi calculating the design error a — a.f 
Then (73) becomes 

„ . . K, + gi2- • • • K^z'" , , 



t The truriciUed series is merely the polynomial on the left side of (56) which 
would lie olitiiiiied if the filter selectivity were ignored, and m were given the 
maximum value, n. 
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If a sharp cut-off were not required, this approximate H{z) could be 
made exactly zero, by using Kk = Kk for all coefficients. Then the actual 
design error would be determined by the approximation inherent in the 
use of (74) in place of (73). For high selectivity, however, X„ should be 
much larger than Kn , as large as possible within assigned limits on 
a — a in the useful range. (It is readily estabhshed that K^z" will deter- 
mine a at asymptotically high frequencies.) The other Kk are then to 
be adjusted so that a — a exhibits the desired "equal ripples." 

The following H(z) has the functional form (74), and also meets con- 
ditions (45) : 

^ , Ki K^ 

„., ^ .. ^ ^' ""^" (75) 

Multiplying Gz'" into the numerator gives a rational fraction which is 
obviously consistent with (74). The coefficients Kk of (74) which corre- 
spond to (75) are : 

if fc = Z, + GiCn., (76) 

In the useful interval [^ K^/z^''] is [S Kkz'"]*. Hence the polynomials 
in z and l/z have identical magnitudes, in the useful interval; and, 
since also \z\ = 1, | H{z) \ = | G | in (75). The phase variation, over 
the useful interval, is the same for H{z) as for z^", which yields the same 
number of ripples in a — o: as an ordinary Tchebycheff filter of like 

degree, t 

The constant G is arbitrary, except that its sign must be properly 
chosen to avoid non-physical natural modes. Increasing G in(;reases the 
filter selectivity, but also increases a — a in the useful interval. G and n 
are to be chosen together, to realize an assigned selectivity within an 
assigned limit on distortion. 

The above analysis may be related to the following filter problem: 
Required to design a filter which has flat gain, in the useful interval, 
but which has m assigned frequencies of infinite loss, in addition to n 
arbitrary natural modes (m ^ n). The n arbitrary natural modes may 
be regarded as compensating for gain variations due to the assigned 
frequencies of infinite loss, in the useful interval, while reinforcing their 
effects at other frequencies. Compensation of effects of the infinite loss 
points is the same as simiilation of the effects of natural modes at the 
same (assigned) frequencies. The approximation in the useful interval 

t This assumes an a with the general characteristics described in Seetioii 8, 
which are such that the numerator and denominator of the fraction in (75) will 
each have a net phase shift of zero, across the useful interval. 
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is to be no better than necessary, so that there may be a maximum 
reinforcing of losses at other frequencies. In these terms, (58), and 
/, KkZ' m (73), correspond to the assigned natural modes (at the 
same locations as the assigned frequencies of infinite loss). Then the 
ideal ^ KkZ''' is itself a polynomial, of degree m ^ 71, and (74) is exact, 
rather than an approximation to (73). Then (76) determines the n 
arbitrary modes in such a way that the net filter gain approximates 
zero in the Tchebycheff sense, over the useful interval. 

A different procedure for obtaining the same result is described in the 
author's paper "Synthesis of lleactanee 4-Poles".^ The above analysis 
of the filter problem is of interest in relating the more general synthesis 
techniques, in terms of Tchebycheff polynomial series, to previous filter 
theory. 

Similar filters have also been obtained by Matthaei^, on a potential 
analogy basis. He includes, however, somewhat more general filter char- 
acteristics, for which he obtains only approximately equal ripple errors. 
Analysis of the sort described above may be used to clarify Matthaei's 
analysis of the conditions under which he obtains exactly equal-ripples. 

Equation (75) may be related to work of Bashkow''. The (arbitrary) 
amplitude of the (equal) maxima of | a — « | , computed from H{z) of 
(75), depends only on [ G [ . The frequencies at which the maxima occur 
correspond to phase H{z) = sr, which is independent of \G \ . Thus, 
the locations of the maxima are invariant to the arbitrary amplitude, 
within the range where (75) applies. (75) applies only when (74) may be 
used in place of (73). Generally, (74) only approximates (73), and the 
approximation introduces small variations in the maxima of | a — a | 
(when a coi-responds to (76)). If the maxima themselves are sufficiently 
small, the small variations will be large percentage variations; and the 
adjustments to compensate for the variations ^\ill yield significant 
shifts in the location of the maxima. In other words, the locations of 
the maxima of | a — a | , required for equal amplitudes, are largely 
invariant to the magnitude of the et^ual amplitudes, but only to an 
approximation which becomes worse as the amplitudes are decreased. 

Baslikow states the invariance of the frequencies of maximum 
\ a — a I , as a more or less empirical conclusion, based on a quite dif- 
ferent approach to the same synthesis problem. 

Equation (75) may be related also to work of Kuh.^ The natural 
modes z<, are zeros of 1 + H(z). In other words, ff (z„) = — 1. Using the 
H{z) of (75) gives the following: 

K, + K,zl + . . ■ Kr^zT = -GzT Ik, + ^ + ■■■—[ (77) 
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It must be remembered, however, that this formulation is permissible 
only if the approximations, inherent in (75), are justified when 2 = z, 
(as well as in the useful interval). Taking the logarithm of each side 
gives: 

log 1^0+ ■■■Kj."} 

( - K^ (78) 

= log i-G) + 2n log z. + logjiiCo + ■ - ■ ^| 

Equation (58) may now be applied, to replace the logarithms by 
power series, provided the truncation of the series is again justified, and 
provided the convergence of X) ^2*2"*' is also proper, at both z = z, and 
z = l/zt, . This gives 

-E C,d' = -E C^^Af + 2n log z. + log i-G) (79) 

(Summations E ^^'^ ^'^ ^^^^^ respect to k; and there is one equation for 
each a.) This is a suitable rule, for obtaining an | a — a | with equal 
maxima, whenever the approximations are in fact luiimportant. It is 
not at all clear, however, just when the approximations become sig- 
nificant. 

Kuh uses the potential analogy approach for the same sort of syn- 
thesis problem. He spaces the natural modes along a p-plane contour 
defined in fairly comphcated potential analogy terms. It can be shown, 
however, that mapping his potentials from p plane to z plane leads to 

(79). 

When the network is to approximate a in the useful interval, but is 
not required to supply selectivity at other frequencies, the approxima- 
tion (74) is usually untenable. It is generally necessary to retain the 
exact formulation (73). 

When selectivity is not required, the phase excursion of H{z), m the 
useful interval, can usually be increased to that of z^"^~ (as in (G7), 
corresponding to Gt = C.t , fc ^ n). As a step toward meeting the first 
condition of (45), one may then write (73) as follows: 

E^ n+l 1 

-- — z + 2^ Qk 2fc 

m^ ^ ,2-1+2 -^"+1 1 ^ 

z) = -Kn+iz J^K^'" (80) 

Q, = -— , fc = 1, ■ ■ ■ , n + 1 

The coefficients Kk and ^i^" are fixed by a. The only arbitrary 
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design constants are the Qk . They are to be small enough so that they 
do not affect the total phase excursion H(z), when \ z \ = 1. Their 
specific values are to make | H{z) \ approximately constant, when 
I 2 I = 1. In general the (required) series in z^ in the numerator makes 
it extremely difficult to determine the rec^uired values for the Qk ■ No 
reasonably simple general procedure has yet been found. 

20. ARBITRARY RATIONAL FRACTIONS 

The preceding sections were devoted to the approximation a to a, 
using n arbitrary natural modes, but no arbitrary frequencies of in- 
finite loss. Similar techniques may be used when there are n" arbitrary 
natural modes and n' arbitrary frequencies of infinite loss. As the appli- 
cations become more involved, however, routine calculations must be 
supplemented increasingly with an element of art. 

For simultaneous design of natural modes and frequencies of infinite 
loss, we must go back from (31) to the a formulation in (21). This we 
shall now write; 

The functions A'' and D are polynomials. The coeflficients will be defined 
as follows: 

A'^ = 7vii + Ai2 ■ • ■ K„,,z 

D = \ + K[z' ■■■ K'^J"' 

By comparison with (21), the zeros of .V, in terms of z^, are the z'J'. 
(Note the minus sign in 81.) The zeros of D are then the zl'. For physical 
networks, n" ^ n'. 

E(iuations (50), de.scrihing a, may be retained as they stand. Com- 
bining (50) and (81) gives, in place of (52): 

a^ a=j: (C, - C,,)T2k 

^ - „ rv 1 (83) 

i: (Cu- - Cn)r'' = -log j^^ R{z)R{-z)^ 

The function R{z)R{ — z) is exactly the same as before. The reciprocal 
of the function will still he ^ K^z^'' , with the Kk related to C-ik by (58), 
(59), (01). The new rational fraction N/D will appear where previously 
we had the polynomial A'n + ■ ■ ■ KnZ^" . 
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Accordingly, the rule for C^h = C-ih ,k^m, now becomes, in place of 

(56), 

~ = J:K,^"' (84) 

This condition may be used to determine the coefficients Kk , Kk of 
A'" and D (in combination with conditions of other sorts, when 
m < n" •{■ n'). When the coefficients have been calculated, the (j-plane) 
natural modes z'J may be determined from the roots of A'', exactly as 
the 2o of previous sections. The infinite loss points z„ may be calculated 
from the roots of D, in exactly the same way except that Re z, need 
not be negative. Signs of the zl must be such that complex and imaginary 
2„ are in conjugate pairs. Note that there can be conjugate imaginary z, 
only if D has a corresponding double negative real zero. 

When m = n" + n', the following method may be used to calculate 
the K'I and K[ determined by (84). The relation is first multiplied by 
/), to get : 

Then algebraic manipulation is used to evaluate the power series equiv- 
alent of the right hand side, through terms of order n" + n', using the 
known values of the Kk , but general values of the K^ of D. Each coef- 
ficient is a linear function of the unknowns, Kk ■ 

Now the polynomial A^, in (85), has no terms of order k > n" . There- 
fore (85) requires zero coefficients in the expansion of the right hand 
side, from order n" + 1 to order n" + n' . Equating these coefficients 
to zero gives n' linear equations in the n' unknown Kk . Solving for the 
Kk determines polynomial D. The values calculated for the Kk may 
then be used in lower order coefficients of the expansion of the right 
hand side of (85), which are exactly the coefficients Kk of A'^. 

When n" — n' = or 1, a continued fraction method is likely to be 
preferable. Various established techniquesf may be used to convert the 
series ^ KkZ^ into a continued fraction of the form: 

E Kkz'' = ao + 



ai . 1 ^ 

? + ~^ i (86) 



a2 + 



- + -V- 



t See, for example, Fry's applications of continued fractions to network de- 
sign. " 
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If the continued fraction is truncated after the term of order m, and is 
rearranged as a rational fraction N/D, it \v\\\ obey equation (84). The 
degrees of N and D will be such that n" + n' = m, and n" — n' = 
or 1. The continued fraction may be associated with the hypothetical 
ladder network shown in Fig. 10, with variable impedance shunt branches 
proportional to z . The impedance of the (truncated) ladder is N/D. 



21. ACCURACY 

The accuracy of match, a — a, may again be evaluated from the 
final network singularities; or by (41), Tvdth H(z) as in (30), before the 
singularities have l)een determined from roots of N and D. A rougher 
estimate may again be obtained from the error in the first immatched 
term in the Tchebycheff polynomial series. As before, (equation (72)), 
this is equal to the leading term in H(z), \\ith z^'"'^~ replaced by — T2m+2 ■ 

The rational fraction in (30) is the same as our present N/D. In terms 
oiN/D, (30) becomes: 



^Riz)R{-z) = l-\- Hiz) 



(87) 



If (86), or Fig. (10), is used to determine A'' and D, the leading term in 
H{z) turns out to l)c: 



H{z) = 



(-)"'+V'"+^ 



(aiOa ■ • • a„,ya,r^iKo 



The corresponding mismatch in Tchebycheff polynomial terms is: 

(-)- 



C2m+2 — C;2m+2 — 



(aia2 ■ ■ • amYa„^iKo 



(88) 



(89) 



22. ZER08 AND POLES 

When frequencies of infinite lo.ss are to be chosen, as well as natural 
modes, the situation in regard to | 2o | < 1 is somewhat less favorable. 




Fig. 10— A ladder network representation of the continued fraction (86). 



650 THE BELL SYSTEM TECHNICAL JOURNAL, JULY 1952 

It is Still true that 1 + H{z) mil have the same number of zeros as 
poles in the region \z\ < 1, so long as a — 5 is reasonably small in the 
interval \z\ = 1. In equation (87), however, the poles of 1 + H{z) 
include the zeros zl of D (the arbitrary infinite loss points), as well as 
the poles of R{z)R{-z). When the coefficients of D are to be chosen as 
in the previous section, the contour rule merely says that any z^ and z^ 
m the region ] 2 [ < 1 will occur in like numbers. 

Fortunately, the frequent occurrence of [ z„ | < 1 is softened by the 
following curious circumstance. Almost always, any | 2„ ] and | z„ | < 1 
are so nearly identical that factors (z — z[) and (z — z") may be can- 
celled out without any important effect on H{z), or a — a. Cancellations 
of this sort were encountered a number of times before an explanation 
was discovered. Actually the explanation is quite simple. 

At any zero of 1 + H{z), H{z) = -I. On the other hand, H{z) is 
small when \z\ = 1. Generally, it is much smaller in most of the inter- 
val I z I < 1. For instance, when Co* = Cat through Ic = m, /f(z) is pro- 
portional to ^^"'■^^ in the neighborhood of z = 0. As a result, 1 H(z) | 
rarely becomes as large as 1, in the region \z\ < 1, except in the very 
close proximity of a pole. In other words, in the region | z | < 1, any 
zero z'J is usually very close to a pole z„ — usually so close that the 
correspoTiding factors z — z„ may be cancleled out without significant 
effect on a — a. 

The occurrence of non-physical natural modes (Rez'J = 0) is the same 
as before; but adjustments to correct for these, in an efficient manner, 
are much more complicated. In addition, there may be non-physical 
infinite loss points, zl . To correct for non-physical singularities, the 
simplest procedure would be to change one or both of the highest order 
coefficients in A^ and D of (82), that is /-:'„',, and /<:'„' . Th|s would be 
inefficient, however, for it would spoil the match of C„" to €„■' , or €„• 
to C„' . The unmodified design, defined by (84), can match terms through 
order n" -\- n', and it is desirable to change oidy the highest order terms 
in adjusting the design. 

More efficient adjustments are in fact feasible. They sometimes re- 
quire an increased element of art; hut the art may be based on specific 
principles. Some particularly useful principles are described in the next 
section. These apply to various other corrections besides correction of 
non-physical zeros and poles. Examples arc reduction in phase to make 
two-terminal realization possible, and increase in shunt capacity in 
two-terminal designs. In general, they offer a means of making 
m < n' -\- n" in (84), and using the remaining degrees of freedom to 
meet other conditions. 
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23. MODIFICATION OF N/D 

Suppose Ni/Di and N^/Di are two rational fractions, representing 
special choices of A'^ and D in (82), such that: 

(90) 

Then suppose F is a fiuiction of z such that 

F = (91) 

The following combination represents an approximation to "^ Kkz'', 
of the order indicated:! 

m = Till , or 7712 + M, whichever is the smaller. 

The left hand side of (92) may be used as N/D in (84), to match 
C-ik to C'ik through /; = m. Adjustment of the function F may be used to 
satisfy other requirements, in addition to accuracy specifications. 

Frequently, Ni/Di may be the rational fraction corresponding to 
truncation of the continued fraction, (86), after the term in an"+n' . 
Then Ni/Di is likely to be a truncation of order m < n" -\- n'. The 
corresponding F is likely to be proportional to z^'', or at most a simple 
polynomial in z'. When F is a constant, and m = n" + n' — 1, the 
use of these particular rational functions in (92), to determine N JD, 
corresponds to matching C«k to C»k through k = n" -\- n' — 1, but 
leaving C'-.(„"+„') subject to adjustment. Specificially, C'i(n--^n-) depends 
on the choice of F, which may hinge upon such special conditions as 
the elimination of non-physical singularities. 

Problems which call for more complicated combinations are by no 
means uncommon. Skill may be needed in the choice of specific com- 
binations which will solve specific problems. Computations may be 



t The relatioDship is easily established by noting that: 

— - 2 A'i-z" — - 2 Kiz^ 

"' + ""' - V AV" = '^ - + F ^ - (93) 

D, + FD. ' D. D, ^■'" 
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systematized, to a considerable extent, by using the error formula (88), 
and other relations between the coefficients of the continued fraction, 
and the rational fraction truncations of various orders, f 

24. APPROXIMATION TO BOTH GAIN AND PHASE 

The applications described in previous sections relate to approxima- 
tions to prescribed gain, a, mthout regard to the associated phase. 
Quite similar methods apply, liowever, to the simultaneous approxima- 
tion of gain and phase. 

The starting point is equation (20). Replacing products of factors by 
polynomials gives, in place of (81) : 



a -|- j/3 = ^ CkTk , k even and odd 



(94) 



The polynomials are now as follows, in place of (82): 

(95) 



■AT T^" 1 Trf T7-" ^' 



D = l-\r K[z--- iC'z"" 

(If only natural modes are to be used, the suitable replacement for the 
first equation or (55) is here obtained merely by using D ^ 1, and K^ , 
z, Za , in place of their squares.) 

A comparable expression is needed for the assigned gain and phase 
a -\- i^. In place of (50), we now repeat (22), and redefine the coeffi- 
cients iCk in accordance with 

a -{- iB = 'y.CkTk , k even and odd 

^ - , (96) 

D hC,z' = log R{z) = -log Z Kkz' 

The definition of Kk has been changed in such a way that it is now 
related to Ca:/2 exactly as it was pre\aously (in 58) related to C^k . 
Equations (59), (61) may be appliedto calculating the Kk by merely 
substituting therein a Ct/2 for every Cat . 



t For example, a simple recursion formula may be used to assemble the poly- 
nomials N and D which correspond to truncation of the continued fraction (86) 

at a number of different points. Specifically, F„ = P„-i + ^-^ Pn-a , where P 

is either N or D and P„ corresponds to truncation of the continued fraction after 
the term in On . The formula holds for i £ 2. 
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Equations (84), (85), (86) may now be modified, for the new A'', D 
and Kk , by merely using z in place of z wherever it occurs (including 
2 in place of 2^*).t The modifications of equations (84), (85), and the 
truncation of (8(i) after a^ now lead to Ck = Ct , k ^ m, instead of the 
previous C^k = C^- . This means that m must be tmce as great to match 
coeflficients out to the same actual orders. This is to be expected since 
now one half of our design parameters are used to approximate phase 
$, leaving only half for approximating gain a. Equation (89) must be 
changed not only in regard to the orders of Ck , Ck , but also in regard 
to the factors ^ in (94), (96). This gives 

Cm+l — Cm+l _ (~) ' /q-x 

2 (aiai ■ ■ ■ a„,y a^+iKo 

The most important change is in regard to the zeros and poles z„ . 
The polynomials jV and D now determine z'J and si directly, instead of 
their squares. There is no opportunity to adjust the sign of Re z'J by 
choosing the correct sign of s/zT'- When non-physical singularities ap- 
pear, adjustments of high order coefficients may be tried. Section 23 
applies provided z'^ is replaced by z. If the specification of the problem 
permits added delay, linear phase may be added to a -(- 'i^ to increase 
the probability of physical singularities! . (Addition of linear phase 
changes only Ci , in ^ CkTk ■ A negalive change in Ci increases the 
delay.) 

2i). APPROXIMATION TO AN ASSIGNED PHASE |9§ 

Sometimes it is required to approximate an assigned phase, without 
regard to gain. More commonly, it is required to approximate an as- 
signed phase, using an "all-pass" network, which has a theoretically 
zero gain. These two problems, however, are very nearly identical, due 
to circumstances explained at the end of this section. 

For approximation to phase only, we go back to the /3 equation in 
(21). As before, products of factors (z — z„) are replaced by polynomial 

m _ me 

t and = in place of = . 

X Tlie well known relation between the gain and phase of any physical network 
(See for instance Bode'") may give some information regarding the reasonable- 
ness of $. It must be remembered, however, that departures of network gain a, 
from the assigned gain a, outside the useful interval, may affect the permissible 
phase 3, within the interval. 

§ Up to the present, applications to phase problems have not been developed 
to the game extent as for gain. Techniques have been explored, however, to deter- 
mine how such applications may in fact be carried out. 
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equivalents. Then, in place of (81) or (94), we have 

1^ = E c,n 



E C^' = -log ^ 



k odd (98) 



(99) 



Using n to represent n" + n', the total number of network singularities, 
we may write A'^ and D as follows, in place of (82) or (95) : 

N = I + KiZ-\- K^i- + iCss' ■ • • + Ky 

D = I - K,z + Kiz - aV ■ ■ ■ i-YK^z" 

Notice that A'' and D are related by 

D{z) = N{-z), (100) 

which is required by the form of the /? equation in (21). 

To arrive at a design procedure most easily (but not the simplest 
design procedure), one may express the assigned gain ^ in the following 
way (comparable to (58) and (96)): 

^^ = E C,T, , fc odd ^^^^^ 

Y.C.z' = -logE^*-^' 

Coefficients Xt- may again be calculated by a modification of (61). This 
time Cat is replaced by Ck , wherever it appears in (61), and thenall 
even ordered Ck are made zero (since only odd terms appears in E ^t^;*" 
of (101)). Note that even ordered K^ are not usually zero, even though 
even ordered Ck are. 

The degrees of A'^ and D, in (99), are such that we can make Ck = Ck 
through terms of order h = 2n. This requires merely: 



2n 



^ = E^^' 002) 

As stated, the condition applies to Ck of both even and odd orders. 
Since even ordered Ck are zero, it means that at least n even ordered 
Ck will be zero, in addition to the match between n odd orders. (102) 
is sufficient to determine an N and a D without reference to (100). If 
the (equal) degrees n of (99) are assumed, however, the A^ and D deter- 
mined by (102) will be found to obey (100) automatically (provided 
E Kkz" corresponds to an odd series E C't^*. ^^ here assumed). f 
A simpler method for computing the same AT and D takes advantage 

t This was discovered by Mrs. M. D. Stoughton. 
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of the known relation (100), connecting N and D. T^et E and be re- 
spectively the sums of even and odd terms in N. Then N in K -\- and 
(100) requires D to be E - 0. The ratio 0/E may be related to 2 ^'^2:* 
of (98) as follows: 

I = -tanh|j:C,3* (103) 

Now let two convergent series, respectively even and odd, be such 

that: 

g = -tanhiD^^" (10-1) 

Let coefficients Kt be defined by: 

E + = E kW (105) 

Then E and are respectively the sums of the even and odd terms. 
The complete series may now be related to the (odd) series ^ Ckz'' as 
follows : 

X hCz = -log i: Klz' (106) 

This fixes the Ri oi E + in terms of the Ck . 

To make Ck = Ck through //; odd orders, (102) is now replaced by 

m 

mo 

The symbol = designates equality of power series through m odd orders. 
(.4/^ even terms are now zero on both sides.) The right hand side may be 
expressed as a continued fraction of the following form, comparable to 

(86): 

E 

(108) 



Truncation after only the m*'' term gives the 0/E of (107). 

The coefficients of and E may be calculated by an appropriate 
modification of (01). (Calculate like Kk of (101), after dividing all Ck 
by 2). After and E have bcett evaluated, by truncating the continued 
fraction (108), their sum gives polynomial N of (99). 
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The natural modes and frequencies of infinite loss are determined 
from the zeros of the polynomial N. By (21), each zero is either a (z- 
plane) natural mode, zj , or the negative of an infinite loss point —z^ . 
If gain variations are inconsequential, there is likely to be some latitude 
in designating each zero as a ^^ , or as a — z„ . 

A zero of A'' with a positive real part would make a non-physical natural 
mode, and hence it must be a —z'^, corresponding to an infinite loss 
point. A zero of A'^ witli a negative real part can be a natural mode Za , 
but this may not be required. It may be either a z„ or a —z, , pro\'ided 
conjugate zeros arc assigned in the same way, and provnded the total 
number of —z'g does not exceed the total number of z'J . The latter con- 
dition requires: 

At least half the zeros of N 
must have negative real parts. 

The continued fraction (108) shows how many zeros will have nega- 
tive real parts, before any zeros have been calciUated. The following 
theorem makes this easy : 

The number of zeros of N which have negative real parts 
is equal to the number of positive coefficients in the trun- 
cation of the continued fraction (108) which determined N. 

When gain is not to be disregarded, but is to be exactly zero, the 
synthesis technique needs few changes. The phase of an "all-pass" net- 
work is related to the natural modes z, as follows: 

i^ = 'ZCkTu, fcodd 






^ ^ = -log ) Zi ^1^) 



This may be regarded as a special case of (20) for a = (wiiich makes 
Ck = for k even, and also happens to require 2, = —z^)- In functional 
form however, it is more like i^ of (21). It differs in only two regards. 
In the power series in z, each Ck is divided by two. In the rational frac- 
tion in z, all the zeros correspond to natural modes, and the poles cor- 
respond to frequencies of infuiite loss ; but the poles are also exactly the 
negatives of the zeros, as in the ip equations of (21). 

Accordingly, the phase synthesis technique which ignores gain varia- 
tions may be applied to the zero gain form of the problem by cutting 
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every Ck in two. All zeros of N = E + must be construed as natural 
modes z, . Finally, the network must have as many infinite loss points 
as natural modes, such that z„ = —z^. (Integer n is now the number 
of natural modes, rather than the total number of singularities.) 

For physical networks, all the first n terms of the continued fraction 
(108) must now be positive. To meet this condition it may be necessary 
to add linear phase to the assigned phase (by adding a negative cor- 
rection to Ci). It appears that sufficient linear phase will always lead 
to a physical design, provided the number of modes n is increased to 
retain a reasonable accuracy. 

26. LINEAR PHASE 

When the assigned phase ^ is linear, the calculations are relatively 
simple. 

If a delay D is to he approximated over a frequency interval extending 

to CO = COc , 

i^ = -D(^,Ti (110) 

If delay D is to be realized without regard to gain variations, the ap- 
propriate 0/E is 

2 = Unh ^ (111) 

A known continued fraction expansion of tanh X may be applied to 
(111), to obtain the eoefficioiits of (108) without bothering with (105). f 
The result may be arranged as follows: 



(112) 



Truncation of the continued fraction gives 0/E, and then + E. The 

zeros Zg turn out to be proportional to yr , and therefore root extraction 

techniques are required only for one D, for each n. The zeros are tabu- 
lated for sample n's, in Table II. 

t For tlic eximnsion of fanh A', reference may be made to a text on continued 
fractions by Wall", page 3-li), equation 91.6. 
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Table II — Z-Plane Natural Modes for Linear Phase 


n 


DucZa 


1 


-2 


2 


-3 ± i VS 


3 


^4.64438 

-3.67782 ± i 3.50876 



4 -5.79242 ± i 1.73446 
-4.20758 ± i 5.31484 

5 -7.29348 

-6.70392 ± i 3.48532 » 
-4.64934 ± I 7.14204 

6 -8.49668 ± i 1.73510 

-7.47142 ± i" 5.25256 
-5.03190 ± i 8.98532 

The error in the first mismatched Tchebycheff coefficient is a rough 
measure of accuracy. It may be shown to be 

^ r - (-)"(Z)co,)'"^' ,..^. 

C,„+i - C2.+1 - 4„(^.3.g ... (2,, _ i)f(2n + 1) ''''''^ 

This measure of accuracy is plotted in Fig. 11, for various numbers of 
natural modes n. A sample detailed curve of /3 — ^ is shown in Fig. 12, 
with dotted lines corresponding to the estimated error (113). 

If delay D is such that the error is reasonable, all the zeros may be 
natural modes. If these are combined with a like number of infinite 
loss points, such that zl = —z'J , an all-pass network will be obtained, 
instead of one which approximates D without regard to gain. The all 
pass network will produce twice the delay, and twice the nonlinearity 
of phase. In other words, for an all pass network, both coordinates in 
Fig. 11 must be doubled. 

27. SIMPLIFICATION OF SINGULARITY ARRAYS 

In complex communication systems, a single equalizer may be re- 
quired to correct for a number of effects. In a coaxial cable system, for 
instance, a single network in the standard repeater may be required to 
compensate for the following: Cable attenuation, characteristics of input 
an(i output networks, effects of interstages (significant because the feed- 
back is limited), and distortion due to variable controls at mean settings. 
Tchebycheff polynomial methods may not be efficient when applied 
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Fig. 11 — Estimated error for n natural niodea approximating linear phase. 
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directly to all these effects. They may still be useful, however, when 
applied in the follomng way: 

Separate arrays of singularities are determined, which match the 
separate effects to required accuracy, using any convenient methods. 
Minimum network complexity is not required at this point. Combining 
all the singularities in a single array gives an initial design which is 
sufficiently accurate, but may use many more singularities than are 
actually necessary. Tchebycheff polynomial methods are now used to 
obtain a simpler set of singularities, which approximates the initial set 
to sufficient accuracy. This has been designated "boiling down" the 
original set. 

In a problem of this sort the assigned characteristic has the network 
form, as well as the network characteristic which is to approximate it. 
(The example discussed in Sections 10 to 14 is also a problem of this 
sort.) As a result, equations (20) and (21) apply to a and j3, as well as 
to a and p. This makes it possible to replace 22 ^kZ^'^ and 2^ ■^t'^ > of 
(55), (56), (96) etc., by a finite rational fraction N/D. If both a and 3 
are to be approximated, the following is derived from (20). 



n(i-| 



N 



2 ICz' = log K. ) f( = -log ^ (114) 

The singularities i'J , zl correspond, of course, to the network singu- 
larities which are to be boiled down. If only a, or only /S, is to be ap- 
proximated, suitable modifications are readily derived from (21). 
The boiling down is accomplished by requiring 

f = 5 ("« 

where N/D is of lower total degree than N/D. If m = n" + n\ and 
n" — n' = or I, the continued fraction method can again be used. 
This requires expansion of N/D in continued fraction form, instead of 

An example of a boiled down set of singularities is illustrated in Fig. 13. 

28. GENERALIZATION OF THE USEFUL INTERVAL 

All the previous analysis applies to a useful frequency interval — w^ < 
01 < -fojc . Its important characteristics are as follows: It is a single 
continuous interval, with cj = at its center. Useful intervals with other 
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other characteristics may be obtained, within Umits, by changing the 
definitions of z and z„ , in terms of p and p* (equations (5) and (8) ) . 

In all cases, the definition of Tchebycheff polynomial Tk remains 
the same in terms of z. The interval of orthogonality remains \z\ = 1; 
and the relation between p and z is always such that the useful fre- 
quency interval is the p-plane mapping of [ s | = 1. The relation must 
also be such that rational functions of p may be expressed as products 
of rational functions, in z and Ijz respectively, corresponding to (13), 
(14). At the same time, the physical restrictions on network singulari- 
ties po must translate into manageable restrictions on the 2-plane singu- 
larities z^ . It is restrictions such as these that limit the manageable 
useful intervals. 

It is easy to apply the "low-pass" techniques to "high-pass" in- 
tervals, extending from w^ , through co , to — oj^ . The appropriate trans- 
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formation, replacing (5), is: 

^"1 (116) 

2 — — 

Z 

With this transformation, the whole of the previous s-plane analysis 
may be applied at once to high-pass useful intervals, except where 
linear phase is involved. 

The situation is more complicated in regard to "band-pass" intervals. 
If the useful interval includes the frequencies between Uci and 0)^2 , the 
complete useful interval (p-plane mapping of | 2 ] =^ 1) must include 
also the "image" frequencies, between —0:^1 and —Uci- Otherwise, 
conjugate complex 2-plane singularities z^, will not lead to conjugate 
network singularities p, . Whcii there are two disjoint parts of the use- 
ful interval, the appropriate relation between p and z is relatively com- 
plicated. Up to the present, no corresponding technique has been dis- 
covered for approximating assigned phases over band-pass intervals, 
in Tchebyeheff polynomial terms. Gain approximations can be handled, 
however, and for a quite simple reason. Gain functions are even func- 
tions, and behave in the p'^ plane much as gain-and-phase functions 
behave in the p plane. In the p^ plane, —w and -\-oi are identical, and a 
band-pass useful interval is a single segment of the a; axis. 

For gain approximations over a band-pass interval, (5) may be re- 
placed by: 

a + h(z - - 

2- \ 2/ (117) 



p = 



1 + c 



(-;)' 



The three coefficients, a, h, c are subject to two conditions, stemming 
from the requirement that the interval \z\ =^ I must map onto the 
interval coL < w^ < i^li . This leaves one arbitrary degree of freedom. 
Its choice may be related to ordhiary least squares approximations in 
the following way : 

If a = XI C^kT2k , the first n terms approximate a in the least squares 
sense. In other words, the integrated square of the error is a minimum, 
relative to all possible choices of the first n coefficients C2k , provided 
the integration extends over the useful frequency interval, and in- 
cludes an appropriate "weight function". When (117) relates z to p, 
the arbitrary degree of freedom in the choice of the constants a, h, c 
permits selection of any one of a family of weight functions. Conventional 
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least squares analysis may be applied to determine these func- 
tions.! 

In applying least squares analysis, it must be borne in mind that the 
iietworlv gain a does not approximate the assigned gain a in the simple 
least squares sense. When C-n- ^ Cok for k ^ n, a — a depends upon 
two least squares approximations. The first n terms of ^ CikTik represent 
a least squares approximation to a, and are made identical with the 
first n terms of ^ Cik-T^k , which represent a least squares approximation 

to a. 

Wlien (117) relates p to z, 2-plane singularities z„ may be defined by: 



(118) 
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An additional singularity, Zo , m also needed, corresponding to the finite 
poles of (117). It may be defined as follows: 

1 + c ( ^0 - - ) =0 

V 20/ (119) 

I 5(i I > 1 

When Pff — V, ''^ "■ of (2), is expressed in terms of z and z^ , (117) 
introduces denominator factors (1 — z' Izt) and (1 — Xjzisz). As a re- 
sult, a. of (21) must be changed to the follomng, for band-pass intervals: 



a = Z] ^ikT% 



n 1- 



When definite values have been chosen for a, b, c of (117) (in order 
that the Ca- may be calculated), (1 — sVso) in (120) is not subject to 
arbitrary adjustment. This situation can be liandled by defining N/D 
as the rational fraction in the a equations of (21), as before, but re- 



t For generftl di3CUBBion.s of orthogonal functions u,nd least squares approxi- 
mations, see Couraiit arid IIill>ert*, and also a sliorl text by Jackson.'^ 
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placing (84) by 



^-~[^~ 4 



T.K^"' 



(121) 



Fig. 14 illustrates an application of the technique to the simulation of 
a coaxial cable attenuation (which is nearly proportional to -s/u). 



29. RECAPITULATION 

Tchebycheff polynomial series may be applied advantageously to a 
very mde range of network synthesis applications. The scope of their 
usefulness may depend upon the skill of the designer, as with any 
synthesis tools, but the underlying principles are reasonably simple. 
The most important principles are perhaps the following : 

A Tchebycheff polynomial series in frequency may be related to a 
power series in a new variable z. When the Tchebycheff polynomial 
series corresponds to a finite network gain or phase, the power series 
corresponds to an analytic function of z, quite similar in form to the 
network function of p, with singularities at 2;-plane mappings of the 
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Fig. 14 — Simulation of a coaxial cable attenuation — Attenuation at top useful 
frequency = 46 db; Network = four constant-resistance sections. 
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network singularities. Tliis makes it possible to apply power series ap- 
proximation methods, in terms of z, to obtain approximations based on 
Tchebycheff polynomial series, in terms of frequency. 

"Maximally flat" approximations in tei'ms of z may be used to match 
the first m terms in the Tchebycheff polynomial series representing net- 
work gain or phase to the corresponding terms in the series representing 
assigned gain or phase. In this way, a Tchebycheff poljoiomial type of 
least squares approximation to the network function is made identical 
to the corresponding least scjuares approximation to the ideal function. 
The overall error, network function minus ideal function, is then the 
difference between the two least squares errors. 

The 3-plane analysis may also be manipulated, in a quite different 
way, to approach an equal ripple type of approximation (which usually 
represents approximation in the Tchebycheff sense). The complications 
are such that applications have been limited to problems of certain 
quite special types. On the otlier hand, analysis of this sort has been 
found useful in clarifj-ing various other ways of seeking equal ripple 
approximations . 
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